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Description 

Field of the Invention 

5 [0001] The invention generally relates to techniques for monitoring the effectiveness of medical therapies and dosage 
formulations, and in particular to techniques for monitoring therapy effectiveness using viral load measurements. 

Background of the Invention 

10 [0002] It is often desirable to determine the effectiveness of therapies, such as those directed against viral infections 
including therapies involving individual drugs, combinations of drugs, or other related therapies. One conventional 
technique for monitoring the effectiveness of a viral infection therapy is to measure and track a viral load associated 
with the viral infection, wherein the viral load is a measurement of a number of copies ofthe virus within a given quantity 
of blood, such per milliliter of blood. The therapy is deemed effective if the viral load is decreased as a result of the 

15 therapy. A determination of whether any particular therapy is effective is helpful in determining the appropnate therapy 
for a particular patient and also for determining whether a particular therapy is effective for an entire class of patients^ 
The latter is typically necessary in order to obtain FDA approval of any new drug or medical device therapy. Viral load 
monitoring is also useful for research purposes such as for assessing the effectiveness ofnew antiviral compounds 
determine, for example, whether it is useful to continue developing particular antiviral compounds or to attempt to gain 

20 FDA market approval. L . . 

[0003] A test to determine the viral load can be done with blood drawn from T-cells or from other standard sources 
The viral load is typically reported either as an absolute number, i.e. , the number of virus particles per milliliter of blood 
or on a logarithmic scale. Likewise, decreases in viral load are reported in absolute numbers, logarithmic scales, or as 
nprcentaoes 

25 [0004] It should be noted though that a viral load captures only a fraction of the total virus in the body of the patient, 
i e it tracks only the quantity of circulating virus. However, viral load is an important clinical marker because the quantity 
of circulating virus is the most important factor in determining disease outcome, as changes in the viral load occur pnor 
to changes in other detectable factors, such as CD4 levels. Indeed, a measurement ofthe viral load is rapidly becoming 
the acceptable method for predicting clinical progression of certain diseases such as HIV. 

30 [0005] Insofar as HIV is concerned, HIV-progression studies have indicated a significant correlation between the nsk 
of acquiring AIDS and an initial HIV baseline viral load level. In addition to predicting the risk of disease progression 
viral load testing is useful in predicting the risk of transmission. In this regard, infected individuals with higher viral load 
are more likely to transmit the virus than others. 

[0006] Currently there are several different systems for monitoring viral load including quantit.ve polymerase chain 

35 reaction (PCR) and nucleic acid hybridization. Herein, the term viral load refers to any virological measurement using 
RNA DNA or p24 antigen in plasma. Note that viral RNA is a more sensitive marker than P 24 antigen. P 24 antigen 
has been shown to be detectable in less than 50% asymptomatic individuals. Moreover, levels ofviral RNA rise and 
fall more rapidly than levels of CD4+ lymphocytes. Hence, changes in infection can be detected more quickly using 
viral load studies based upon viral RNA than using CD4 studies. Moreover, viral load values have to date proven to 

40 be an earlier and better predictor of long term patient outcome than CD4-cell counts. Thus, viral load determinations 
are rapidly becoming an important decision aid for anti-retro viral therapy and disease management. Viral load studies, 
however have not yet completely replaced CD+ analysis in part because viral load only monitors the progress of the 
virus during infection whereas CD4+ analysis monitors the immune system directly. Nevertheless, even where CD4+ 
analysis is effective, viral load measurements can supplement information provided by the CD4 counts. For example, 

45 an individual undergoing long term treatment may appear stable based upon the observation of clinical parameters 
and CD4 counts. However, the viral load of the patient may nevertheless be increasing. Hence, a measurement of the 
viral load can potentially assist a physician in determining whether to change therapy despite the appearance of long 
term stability based upon CD4 counts. . 
[0007] Thus, viral load measurements are very useful. However, there remains considerable room for improvement. 

so One problem with current viral load measurements is that the threshold level for detection, i.e., the nadir of detection, 

is about 400-500 copies per m iter. Hence, currently, if the viral load is below 400-500 copies per milliliter, the virus 

is undetectable. The virus may nevertheless remain within the body. Indeed, considerable quantities ofthe virus may 
remain within the lymph system. Accordingly, it would be desirable to provide an improved method for measuring viral 
load which permits viral load levels of less than 400-500 copies per milliliter to be reliably detected. 

55 [0008] Another problem with current viral load measurement techniques is that the techniques are typically only 
effective for detecting exponential changes in viral loads. In other words, current techniques will only reliably detect 
circumstances wherein the viral load increases or decreases by an order of magnitude, such by a factor of 10 In other 
cases viral load measurements only detect a difference between undetectable levels of the virus and detectable levels 
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of the virus As can be appreciated. It would be highly desirable to provide an improved method for tracking changes 
in viral load which does not require an exponential change in the viral load fordetection or which does not require a 
change from an undetectable level to a detectable level. Indeed, with current techniques, an exponents! or sub-expo- 
nential change in the viral load results only in a linear change in the parameters used to measure the viral load. It would 
5 instead be highly desirable to provide a method for monitoring the viral load which converts a linear change in the viral 
load into an exponential change within the parameters being measured to thereby permit very slight vanations in viral 
load to be reliably detected. In other words, current viral load detection techniques are useful only as a qualitative 
estimator, rather than as a quantative estimator. 

[0009] One reason that current viral load measurements do not reliably track small scale fluctuations in the actual 

10 number of viruses is that a significant uncertainty in the measurements often occurs. As a result, individual viral load 
measurements have little statistical significance and a relatively large number of measurements must be made before 
any statistically significant conclusions can be drawn. As can be appreciated it would be desirable to provide a viral 
load detection technique which can reliably measure the viral load such that the statistical error associated with a single 
viral load measurement is relatively low to permit individual viral load measurements to be more effectively exploited. 

15 [0010] Moreover, because individual viral load measurements are not particularly significant when using current 
methods treatment decisions for individual patients based upon the viral load measurements must be based only upon 
long term changes or trends in the viral load resulting in a delay in any decision to change therapy. It would be highly 
desirable to provide an improved method for measuring and tracking viral load such that treatment decisions can be 
made much more quickly based upon short term trends of measured viral loads. 

20 [0011] As noted above, the current nadir of viral load detectability is at 400-500 copies of the virus per milliliter. 
Anything below that level is deemed to be undetectable. Currently the most successful and potent multi-drug therapies 
are able to suppress viral load below that level of detection in about 80-90 percent of patients. Thereafter, viral load is 
no longer an effective indicator of therapy. By providing a viral load monitoring technique which reduces the nadir of 
detectability significantly, the relative effectiveness of different multi-drug therapies can be more effectively compared. 

25 indeed new FDA guidelines for providing accelerated approval of a new drug containing regimen requires that the 
regimen suppress the viral load below the current nadir of detection in about 80 to 85 percent of cases. If the new 
regimen suppresses the viral load to undetectable levels in less than 80 to 85 percent of the cases, the new drug will 
gain accelerated approval only if it has other redeeming qualities such as a preferable dosing regimen (such as only 
once or twice per day), a favorable side effect profile, or a favorable resistance or cross-resistance profile. Thus, the 

30 ability of a regimen to suppress the viral load below the level ofdetection is an important factor in FDA approval. How- 
ever because the level of detectability remains relatively high, full approval is currently not granted by the FDA solely 
based upon the ability of the regimen to suppress the viral load below the minimum level of detection. Rather, for full 
approval the FDA may require a further demonstration of the durability of the regimen, i.e., a demonstration that the 
drug regimen suppresses the viral load below the level of detectability and keeps it below the level of detectability for 
35 some period of time. 

[0012] As can be appreciated, if a new viral load measurement and tracking technique were developed which could 
reliably detect viral load at levels much lower than the current nadir of detectability, the FDA may be able, using the 
new technique, to much more precisely determine the effectiveness of a drug regimen for the purposes of granting 
approval such that a demonstration of the redeeming evalities will no longer be necessary. 
40 [0013] For all of these reasons, it would be highly desirable to provide an improved technique for measuring and 
tracking viral load capable of providing much more precise and reliable estimates of the viral load and in particular 
capable of reducing the nadir of detectability significantly. The present invention is directed to this end. 

Summary of the Invention 

45 

[0014] In accordance with the invention, a method according to claim 1 is provided for determining the effectiveness 
of a therapy, such as an anti-viral therapy, by analyzing microarray output patterns generated from biological samples 
taken at different sampling times from a patient undergoing the therapy. In accordance with the method, a viral diffusion 
curve associated with a therapy of interest is generated and each of the output patterns representative of hybndization 
so activity is then mapped to coordinates on the viral diffusion curve using fractal filtering. A degree of convergence 
between the mapped coordinates on the viral diffusion curve is determined. Then, a determination is made as to whether 
the therapy of interest has been effective based upon the degree of convergence. 

[0015] In an exemplary embodiment, the viral diffusion curve is spatially parameterized such that samples map to 
coordinates near the curve maxima, if the viral load is increasing (i.e., therapy or dosage is ineffective). In this manner, 
55 any correlation between rate and extent of convergence across different patient samples is exploited to provide a 
quantitative and qualitative estimate of therapy effectiveness. 

[0016] Also in the exemplary embodiment, the biological sample is a DNA sample. The output pattern of the micro- 
array is quantized as a dot spectrogram. The viral diffusion curve is generated by inputting parameters representative 
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of viral load studies for the therapy of interest, generating a preliminary viral diffusion curve based upon the viral load 
studies; and then calibrating a degree of directional causality in the preliminary viral diffusion curve to yield the viral 
diffusion curve. The parameters representative ofthe viral load studies include one or more of baseline viral load (BVL) 
set point measurements at which detection is achieved, BVL at which therapy is recommended and viral load markers 

5 at which dosage therapy is recommended. The step of generating the preliminary viral diffusion curve is performed by 
selecting a canonical equation representative ofthe viral diffusion curve, determining expectation and mean response 
parameters for use in parameterizing the equation selected to represent the viral diffusion curve and parameterizing 
the equation selected to represent the virat diffusion curve to yield the preliminary virat diffusion curve. 
[0017] Also, in the exemplary embodiment, each dot spectrogram is mapped to the viral diffusion curve using fractal 

10 filtering by generating a partitioned iterated fractal system IFS model representative ofthe dot spectrogram, determining 
affine parameters for IPS model, and then mapping the dot spectrogram onto the viral diffusion curve using the IFS. 
Before the dot spectrograms is mapped to the viral diffusion curve, the dot spectrograms are interferometrically en- 
hanced. After the mapping, any uncertainty in the mapped coordinates is compensated for using non-linear information 
filtering. 

15 [001 8] Apparatus embodiment according to claim 1 6 is also provided. By exploiting aspects of the invention, disease 
management decisions related to disease progression, therapy and dosage effectiveness may be made by tracking 
the coordinates on the viral diffusion curve as successive D N A-/RN A-based microarray samples are collected and 
analyzed. 

20 Brief Description of the Drawings 

[0019] The features, objects, and advantages of the present invention will become more apparent from the detailed 
description set forth below when taken in conjunction with the drawings in which like reference characters identify 
correspondingly throughout and wherein: 

25 

FIG. 1 is a flow chart illustrating an exemplary method for determining the effectiveness of a viral therapy in ac- 
cordance with the invention. 

FIG. 2 is a flow chart illustrating an exemplary method for generating Viral Diffusion Curves for use with the method 
of Fig. 1. 

30 FIG. 3 is a flow chart illustrating an exemplary method for mapping dot spectrograms onto Viral Diffusion Curves 

using fractal filtering for use with the method of Fig. 1. 
FIG. 4 is a block diagram illustrating the effect of the fractal filtering of Fig. 3. 

FIG. 5 is a flow chart illustrating an exemplary method for compensating for uncertainty using non-linear information 
filtering for use with the method of Fig. 1. 

35 

Detailed Description of Exemplary Embodiments 

[0020] With reference to the figures, exemplary method embodiments of invention will now be described. The ex- 
emplary method will be described primarily with respect to the determination of changes in viral loads based upon the 
40 output patterns of a hybridized biochip microarray using DNA samples, but principles of the invention may be applied 
to other protein-based samples or to other types of output patterns as well. 

[0021] With reference to FIG. 1, steps will be described for generating viral diffusion curves for use is processing 
DNA biomicroarray output patterns to determine the effectiveness of therapies imposed upon a patient providing sam- 
ples for which the outputs are generated. Then, steps will be described for processing the specific output patterns 
45 using the VDCs. 

[0022] An underlying clinical hypothesis of the exemplary method is that antiviral treatment should inhibit viral rep- 
lication and lower an individual's viral load from baseline or suppress rising values. A stationary or rising viral load after 
the introduction of antiviral therapy indicates a lack of response to the drug(s) or the development of drug resistance. 
The VDC exploits the underlying hypothesis in part by correlating the rate of disease progression to a sample point 
50 value such that a change in sample point indicates progression. 

The Method 

[0023] At step 1 00, parameters representative of viral load studies for the therapy of interest are input. A preliminary 
55 viral diffusion curve is generated, at step 102, based upon the viral load studies. The parameters representative of the 
viral load studies include baseline viral load (BVL) set point measurements at which detection is achieved, BVL at 
which therapy is recommended and viral load markers at which dosage therapy is recommended. At step 104, a degree 
of directional causality in the preliminary viral diffusion curve is calibrated to yield the final viral diffusion curve. 
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[0024] Steps 100-104 are performed off-line for setting up the VDCs. These steps need be performed only once 
for a given therapy and for a given set of baseline viral load measurements. Thereafter, any number of DNA biomicro- 
array output patterns may be processed using the VDCs to determine the effectiveness of the therapy. Preferably, 
VDCs are generated for an entire set of therapies that may be of interest such that, for any new DNA biomicroarray 
5 output pattern, the effectiveness of any of the therapies can be quickly determined using the set of VDCs. In general, 
the aforementioned steps need be repeated only to update the VDCs based upon new and different baseline viral load 
studies or if new therapies of interest need to be considered. 

[0025] In the following, steps will be summarized for processing DNA biomicroarray output patterns using the VDCs 
to determine whether any therapies of interest represented by the VDCs are effective. To determine the effectiveness 
10 of therapy at least two samples of DNA to be analyzed are collected from a patient, preferably taken some time apart, 
and biomicroarray patterns are generated therefrom. In other cases though, the different samples are collected from 
different patients. 

[0026] The output patterns for the DNA biomicroarray are referred to herein as dot spectrograms. A dot spectrogram 
is generated using a DNA biomicroarray for each sample from an N by M DNA biomicroarray. An element of the array 
15 is an "oxel": o(ij). An element of the dot spectrogram is a hixel : h(ij). The dot spectrogram is represented by cell 
amplitudes given by <I>(i,j) for i: 1 to N, and j: 1 to M. 

[0027] Dot spectrograms are generated from the samples taken at different times using a prefabricated DNA biom- 
icroarray at step 106. The dot spectrograms are interferometrically enhanced at step 108. Each dot spectrogram is 
then mapped to coordinates on the viral diffusion curves using fractal filtering at step 110. After the mapping, any 

20 uncertainty in the mapped coordinates is compensated for at step 112 using non-linear information filtering. 

[0028] VDC coordinates are initialized at step 114, then updated in accordance with filtered dot spectrograms at step 
116. A degree of convergence between the mapped coordinates on the viral diffusion curves is then determined at 
step 118 and a determination is made as to whether the therapy of interest has been effective. The determination is 
based upon whether the degree of convergence increases from one DNA sample to another. An increase in degree 

25 of convergence is representative of a lack of effectiveness of the therapy of interest. Hence, if the degree of convergence 
decreases, then execution proceeds to step 120, wherein a signal is output indicating that the therapy is effective. If 
the degree of convergence increases, then execution proceeds to step 122 wherein VDC temporal scale matching is 
performed. Then a determination is made at step 124 whether an effectiveness time scale has been exceeded. If 
exceeded, then a conclusion is drawn that the effectiveness of the viral therapy cannot be established even if more 

30 samples are analyzed. If not exceeded, then execution returns to step 106 wherein another sample taken from the 
same patient at a latter time is analyzed by repeating steps 106 through 118. 

Viral Load Studies 

35 [0029] Viral load studies for therapies of interest are parameterized at step 100 as follows. The therapy of interest 
is selected from a predetermined list of therapies for which viral load studies have been performed. Measurements 
from viral load studies are input for therapy of interest. As noted, the viral load measurements include one or more of 
Baseline Viral Load (BVL) set point measurements at which detection is achieved; BVL at which therapy is recom- 
mended; and VL markers at which dosage change recommended. 

40 [0030] Data for the viral load measurements are obtained, for example, from drug qualification studies on a minimum 
include dosages, viral limits as well as time cycles within which an anti-viral drug is deemed effective. The data is 
typically qualified with age/weight outliers and patient history. Attribute relevant to this claim is the y n or BVL LOW which 
corresponds to the lowest detection limit shown for a therapy to be effective using conventional assays or any other 
diagnostic means. BVL NP . LOW denotes the lowest threshold at which viral load is achieved using a nucleotide probe. 

45 Using interferometric enhancement technique, BVL NP ^ 0W «BVL LOW . 

Generation of Viral Diffusion Curves 

[0031] Referring now to FIG. 2, the viral diffusion curves are generated as follows. An equation is selected for rep- 
50 resenting the VDC at step 200. Expectation and mean response parameters are determined at step 202 for use in 
parameterizing the selected equation. Then the equation selected to represent the VDC is parameterized at step 204 
to yield a numerical representation of the VDC. These steps will now be described in greater detail. 
[0032] These steps populate a canonical machine representation, denoted as VDC(i, V, y, k) (which is a special case 
of Fokker-Planck equation) to calibrate responses from a viral load detection DNA-array based hybridization biomicro- 
55 array. 

i is the index for a diagnostic condition/therapy of interest, 
T denotes the parameter vector characterizing the VDC, 
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X denotes the clinical endpoints vector that indicates detectability thresholds for a specific DNA-hybridization array 
implementation 

k correspond to the uncertainty interval estimates. 
5 [0033] An example of an equation selected for representing the VDC is: 

|e = div( W(x)p) + ^ , p(x.O) = p°(x) 

10 [0034] The potential (x):<R" -> [0, *>) is a smooth function, 0 > 0 is selected constant, and p° (x) is a probability 
density on SR". 

[0035] Preferably, the diffusion potential of the equation and the BVL data are such that: 
15 < c - [BVL NP LOW |BVL LOW ] 

the constant c is generally set to 

[number of amplitude discretization levels] 

20 °~ [log(PCR amplification factor) * avg(oligonucleotides/oxel) * r tagging efficiency! * r binding efficiency - !] 

[0036] Binding efficiency is difficult to quantify analytically for a biomicroarray device technology. Hence, for use in 
the above equation, an estimate of the binding efficiency is preferably employed. A binding efficiency of 30% (.3) is 
appropriate, though other values may alternatively used. Depending upon the specific biomicroarray used, the constant 
25 c typically ranges between 0.0001 to 0.5. 

Expectation and Mean Response Parameters 

[0037] The expectation (u.) and mean response parameters are then determined at step 202 for use in parameterizing 
30 the equation selected to represent the VDC. The expectation and mean response values are determined by: 1 ) per- 
forming conventional PCR amplification; 2) obtaining calibrated viral counts from the PCR amplification; 3) determining 
enhanced and normalized hybridization amplitude mean and variance values corresponding to the calibrated viral 
counts; and 4) matching the enhanced and normalized hybridization amplitude mean and variance values. 
[0038] Two synthetic amplification techniques (in addition to PCR and any designer tagging) are used to achieve VL 
35 estimation above the BVL limit set for the exemplary embodiment of the method, namely (a) readout pre-conditioning, 
and (b) nonlinear interferometric enhancement. Moreover, the expectation match condition implies that: 

Expectation [Log( L interferometrically enhanced imageJ)] v[expressionsetof interest] 
40 Expectation [preconditioned image amplitude] v(expression Mt ^ interest 

= 1 

[0039] Variance matching is done similarly with respect to biomicroarray readout. The lower bound of mean response 
45 value can be given by: 

Variance [Log( L interferometrically enhanced imageJ)] v[expressionsetof interes1] _ 
Variance [preconditioned image amplitude] v[expressjon ^ ^ interest 

50 

[0040] Using the above expression, a conservative lower bound for interferometric enhancement is estimated for 
each nucleotide expression of interest. Since the array fabrication device is assumed to have (i) an identical oligonu- 
cleotide density per oxel and (ii) equal length oligonucleotides, the same mean response amplitude can be assumed. 
If these two assumptions are not met then bounds need to be individually calculated and averaged using the above 
55 formula. Another assumption is that the binding efficiency is statistically independent of the actual oligonucleotide 
sequence. Ifthis assumption does not hold for the specific device technology then the binding efficiency should be 
provided as well for each expressed sequence of interest. So the computational analysis method uses the analytically 
derived lower bounds, as computed using the above equation. This is a one-time calculation only and is done offline 
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at design time. 
Parameterization of the VDC 

[0041] The equation selected to represent the VDC is then parameterized using the expectation and mean response 
values at step 204 to yield a numerical representation of the VDC using 



with constants a, p, y and e « 1 . 
[0042] This utility of this parameterization is established as follows. The VDC canonical representation is based on 
a variational formulation of the Fokker-Planck. The Fokker-Planck (FP) equation, or forward Kolmogorov equation, 
describes the evolution of the probability density for a stochastic process associated with an Ito stochastic differential 
equation. The exemplary method exploits the VDC to model a physical time-dependent phenomena in which random- 
ness plays a major role. The specific variant used herein is one for which the drift term is given by the gradient of a 
potential. For a broad class of potentials (that correspond to statistical variability in therapy response), a time discrete, 
iterative variational is constructed whose solutions converge to the solution of the Fokker-Planck equation. The time- 
step is governed by the Wasserstein metric on probability measures. In this formulation the dynamics may be regarded 
as a gradient flux, or a steepest descent for the free energy with respect to the Wasserstein metric. This parameterization 
draws from theory of stochastic differential equations: wherein a (normalized) solution to a given Fokker-Planck equa- 
tion represents the probability density for the position (or velocity) of a particle whose motion is described by a corre- 
sponding Ito stochastic differential equation (or Langevin equation). The drift coefficient is a gradient. The method 
exploits "designer conditions" on the drift and diffusion coefficients so that the stationary solution of a Fokker-Planck 
equation satisfies a variational principle. It minimizes a certain convex free energy functional over an appropriate ad- 
missible class of probability densities. 

[0043] A physical analogy is to an optimal control problem which is related to the heating of a probe in a kiln. The 
goal is to control the heating process in such a way that the temperature inside the probe follows a certain desired 
temperature profile. The biomolecular analogy is to seek a certain property in the parameterized VDC-namely, an 
exponential jump in the VDC coordinate position for "small linear changes in the viral count". 

[0044] This method is in contradistinction to conventional calibration strategies which obtain a linear or superlinear 
shift in quantization parameter for an exponential shift in actual viral count. 
[0045] As noted, the form of FP equation chosen is : 



where the potential (x): <H" -» [0, «») is a smooth function, p > 0 is selected constant, and p° (x) is a probability 
density on 9?\ The solution p(t,x) is a probability density on 9?" for almost every fixed time t. So the distribution p(t, 
x) £ 0 for almost every (t,x) e (O,«0 X <H°, and 



.for almost every t, e (0,°o). 

[0046] It is reasonably assumed that hybridization array device physics for the DNA biomicroarray (i.e., corresponding 
to the potential y) has an approximately linear response to the nucleotide concentration and the response is monotonic 
with bounded drift. So, 



y(x)=p 0 + p lX + (3 2 x 2 



subject to constraints 




|C = div( W(x)p) + £fi , p(x,0) = p°(x) 



^p(/,x)A=l 



where the partition function Z is given by 
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2 = \eW»dx 



10 



[0047] In this model the basis for device physics design is that the potential needs to be modulated such that it grows 
rapidly enough for Z to be finite. This is not achieved by conventional methods. However, a technique which does 
achieve this result is described in co-pending U.S. Patent Application xx/xxx,xxx, filed contemporaneously herewith, 
15 entitled " Method and Apparatus for Analyzing Hybridized DNA Microarray Patterns Using Resonant Interactions Em- 
ploying Quantum Expressor Functions", which is incorporated by reference herein. 

[0048] The probability measure p s (x)dx is the unique invariant measure for the Markov random field (MRF) fit to the 
empirical viral load data. 

[0049] The method exploits a special dynamical effect to design p. The method restricts the FP equation form above 
20 to a more restricted case: random walk emulating between critical equilibrium points. 

[0050] To aid in understanding this aspect of the invention, consider the diffusion form 

M*-<) = 1 o 2 a2 P(*-0 
3 < 2 a 2 / 

where D 2 = na 2 
and a is constant. 
[0051] A specific VDC shape is parameterized by: 

30 

y(x)=P 0 + Pi x + P2 x2 



subject to constraints 

35 



: = y sin (fi 0 + fi % x + Ax 2 ) j sin cot +*(/) 



with constants a, 0, y and e « 1 . 
[0052] These constants are set based upon the dynamic range expected for the viral load. Thus, if the viral load is 
expected to vary only within a factor of 10, the constants are set accordingly. If the viral load is expected to vary within 
45 a greater range, different constants are employed. The actual values of the constants also depend upon the particular 



50 



55 



[0053] Where the following conditions are met 

Expectation match: E(x) = £^xf(x)dx = /i 

Variance: a 7 = |(x-/j) f t dx 
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And 



]/<*) = 1 



[0054] The expectation and mean response parameters for use in these equations are derived, as described above, 
10 from matching the enhanced and normalized hybridization amplitude mean and variance that correspond to calibrated 
viral counts (via classical PGR amplification). 

[0055] A distribution represented by the above-equations then satisfies the following form with a prescribed proba- 
bility distribution 



x = r sin y(x) I sin cot (/) 



[0056] Assuming: 



and 



y'=^>P>0, p = constanr 



e(0 = e 0 ay 

such that 

a = a 1/3 (y-1)(y + 1)-e 0 a 
and the distribution controlling equation is 

f(x) = 0.5|/(x)| 

such that y(-1)<X<y(1). 
[0057] The characteristic timescale of response for this system is given by 



50 the successive points must show a motion with characteristic timescale. The VDC is designed such that sampling 

time falls well within characteristic time. 

[0058] As noted the actual information used for populating above the parameters is available from the following: 
Baseline viral load (BVL) set point measurements at which detection is achieved; BVL at which therapy is recommend- 
ed; and VL markers at which dosage change recommended. 
55 [0059] The following provides an example of preclinical data that is available to to assist in parameterizing the VDC. 
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NOTIONAL VIRAL LOAD MANAGEMENT EXAMPLE 

[0060] This is a synthetic example to illustrate how data from clinical studies may be used to calibrate the VDC. 
[0061] Viral Load analysis studies, using conventional assays, in HIV progression have shown that neither gender, 
5 age, HCV co-infection, past history of symptomatic HIV-1 infection, duration of HIV-1 infection nor risk group are as- 
sociated with a higher risk of increasing baseline viral load (BVL) to the virologic end-point. However, patients with a 
high (BVL) between 400-6000 copies had a 10-fold higher risk of increasing the level of viral load than patients with a 
BVL below 1500 copies/ml. Thus, baseline viral load set point measurements provide an important indicator for onset 
of disease. 

10 [0062] Initiation of antiretroviral therapy is generally recommended when the CD4 + T-cell count is < 600 cells/mL 
and the viral load level is > 6,000 copies/mL. When the viral load is > 28000 copies/mL, initiation of therapy is recom- 
mended regardless of other laboratory markers and clinical status. 

[0063] Effective antiretroviral treatment may be measured by changes in plasma HIV RNA levels. The ideal end point 
for effective antiviral therapy is to achieve undetectable levels of virus (< 400 copies/mL). A decrease in HIV RNA 
15 levels of at least 0.5 log suggests effective treatment, while a return to pretreatment values (±0.5 log) suggests failure 
of drug treatment. 

[0064] When HIV RNA levels decline initially but return to pretreatment levels, the loss of therapy effectiveness has 
been associated with the presence of drug-resistant HIV strains. 

[0065] The therapy-specific preclinical viral load markers (such as low and high limits in the above example) are 
20 used to establish actual BVL boundaries for the VDC associated with a particular therapy. In this regard, the deterim- 

ination of the BVL parameters is disease specific. For example, in HIV methods such as RT-PCR, bDNA or NASBA 

are used. Other diseases use other assays. Typically, once the parameters of the VDC equation have been set (i.e. 

constants a, p, y and e), only two viral load markers are needed to complete the parameterization ofthe VDC. This is 

in contrast to previous techniques whereby expensive and laborious techiques are required to determine the shape of 
25 a viral diffusion curve. The present invention succeeds in using only two viral load markers in most cases by exploiting 

the canonical VDC described above which has predefined properties and which is predetermined based upon the 

particular biochip being used. 

Calibration of Viral Diffusion Curves 

30 

[0066] Referring again to FIG. 1, the directional causality of the VDC is calibrated at step 104 in the context of an 
NIF discussed in greater detail below. At least arbitrarily selected three sample points are used execute the NIF cali- 
bration computation. The resulting polynomial is used to extracting qualitative coherence properties of the system. 
[0067] The spectral {O] and temporal coherence [ti] is incrementally estimated and computed for each mutation/ 
35 oligonucleotide of interest by a NIF forward estimation computation (described further below). The two estimates are 
normalized and convolved to yield cross-correlation function over time. The shape index (i.e., curvature) of the minima 
is used as a measure for directional causality. Absence of curvature divergence is used to detect high directional 
causality in the system. 

40 Sample Collection 

[0068] Samples are collected at step 106 subject to a sample point collection separation amount. The separation 
amount for two samples is preferably within half a "drug effectiveness mean time" covering 2a population level wherein 
o denotes the standard deviation in period before which a effectiveness for a particular drug is indicated. The following 
45 are some general guidelines for sample preparation for use with the exemplary method: 

It is important that assay specimen requirements be strictly followed to avoid degradation of viral RNA. 

A baseline should be established for each patient with two specimens drawn two to four weeks apart. 

Patients should be monitored periodically, every three to four months or more frequently if therapy is changed. A 

50 viral load level that remains at baseline or a rising level indicates a need for change in therapy. 

Too much significance should not be given to any one viral load result. Only sustained increases or decreases of 
0.001- 0.01 log [conventional methods typically require a 0.3-0.5 log change] or more should be considered sig- 
nificant. Biological and technical variation of up to 0.01 log [typical conventional limit: 0.3 log] is possible. Also, 
recent immunization, opportunistic infections and other conditions may cause transient increases in viral load 

55 levels. 

A new baseline for each patient should be established when changing laboratories or methods. 
[0069] Recommendations for frequency of testing are as follows: 



10 



EP 1 145 177 B1 



establish baseline: 2 measurements, 2 to 4 weeks apart 
every 3 to 4 months or in conjunction with CD4+ T-cell counts 
3 to 4 weeks after initiating or changing antiretroviral therapy 
shorter intervals as critical decisions are made. 

measurements 2-3 weeks apart to determine a baseline measurement. 

repeat every 3-6 months thereafter in conjunction with CD4 counts to monitor viral load and T-cell count, 
avoid viral load measurements for 3-4 weeks following an immunization or within one month of an infection, 
a new baseline for each patient should be established when changing laboratories or methods. 

[0070] The samples are applied to a prefabricated DNA biomicroarray to generate one or more dot spectrograms 
each denoted O(i.j) for i: 1 to N, and j: 1 to M. The first sample is referred to herein as the k=1 sample, the second as 
the k=2 sample, and so on. 

Interferometric Enhancement of the Dot Spectrogram 

[0071] Each dot spectrogram provided by the DNA biomicroarray is filtered at step 108 to yield enhanced dot spec- 
trograms <&(k) either by performing a conventional Nucleic Assay Amplification or by applying preconditioning and 
normalization steps as described in the co-pending Patent Application having Serial Number 09/253,789, entitled 
"Method And Apparatus For Analyzing Hybridized Biochip Patterns Using Resonance Interactions Employing Quantum 
Expressor Functions". The application is incorporated by reference herein, particularly insofar as the descriptions of 
the use of preconditioning and normalization curves are concerned. 

Fractal Filtering 

[0072] Each enhanced dot spectrogram is then mapped to the VDC using fractal filtering at step 110 as shown in 
FIG. 3 by generating a partitioned iterated fractal system 302, determining affine parameters for the IFS 304 and then 
mapping the enhanced dot spectrogram onto the VDC using the IFS, step 306. 
[0073] The VDC representation models a stochastic process given by 




if(x*y) ^ , for some 1 < i£ m; 



»XO(*,>0 = 



otherwise', 



for any (x.y) G<R J and f e p(SR J ) 
[0074] An exemplary partitioned iterated fractal system (IFS) model for the system is 



W«{*, = (|i |l T f »/=1 v 2 m 



where the affine parameters for the IFS transformation are given by 




where the D - origin is given by O^.y^). 
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the R -origin is given by (x^, y' R ) 
spatial transformation matnx is given by o, 
the intensity tilting vector is given by t # - 
the contrast adjuctment is given by y h 
5 the brightness adjustment is given by p /t 

and wherein O represents the enhanced dot spectrogram and wherein u, represents the calculated expectation 
match values 

[0075] This IFS model maps the dot spectrogram to a point on the VDC wherein each VDC coordinate is denoted 
by VDC(t,0) such that 

10 

W[<D, k] -> VDC(k, 0) 

Wherein k represents a sample. 
15 [0076] In the above equation, 0 represent the parameters of the IFS map. 

[0077] Thus the output of step 1 1 0, is a set of VDC coordinates, identified as VDC(k, 0), with one set of coordinates 
for each enhanced dot spectrogram k = 1 , 2...n. 

[0078] The effect of the steps of FIG. 3 is illustrated in FIG. 4 which shows a set of dot spectrograms 450, 451 and 
452 and a VDC 454. As illustrated, each dot spectrogram is mapped to a point on the VDC. Convergence toward a 
20 single point on the VDC implies ineffectiveness of the viral therapy. A convergence test is described below. 

Uncertainty Compensation 

[0079] With reference to FIG. 5, any uncertainty in the coordinates VDC(k, 0) is compensated using Non-linear 
25 Information Filtering as follows. Biomicroarray dispersion coefficients, hybridization process variability values and em- 
pirical variance are determined at step 402. The biomicroarray dispersion coefficients, hybridization process variability 
values and empirical variance are then converted at step 404 to parameters for use in NIF. The NIF is then applied at 
step 406 to the VDC coordinates generated at step 106 of FIG. 1. 

[0080] Nonlinear information filter (NIF), is a nonlinear variant of the Extended Kalman Fitter. A nonlinear system is 
30 considered. Linearizing the state and observation equations, a linear estimator which keeps track oftotal state estimates 
is provided. The linearized parameters and filter equations are expressed in information space. This gives a filter that 
predicts and estimates information about nonlinear state parameters given nonlinear observations and nonlinear sys- 
tem dynamics. 

[0081] The information Filter (IF) is essentially a Kalman Filter expressed in terms of measures of the amount of 
35 about the parameter of interest instead of tracking the states themselves, i.e., track the inverse covariance form of the 
Kalman filter. Information here is in the Fisher sense, i.e. a measure of the amount of information about a parameter 
present in the observations. 

[0082] Uncertainty bars are estimated using NIF algorithm. The parameters depend on biomicroarray dispersion 
coefficients, hybridization process variability and empirical variance indicated in the trial studies. 
40 [0083] One particular advantage of the method of the invention is that it can also be used to capture the dispersion 
from individual to individual, therapy to therapy etc. It is extremely useful and enabling to the method in that it can be 
apriori analytically set to a prechosen value and can be used to control the quality of biomicroarray output mapping to 
VDC coordinates. 

[0084] The biomicroarray dispersion coefficients, hybridization process variability values and empirical variance are 
^5 determined as follows. Palm generator functions are used to capture stochastic variability in hybridization binding 
efficacy. This method draws upon results in stochastic integral geometry and geometric probability theory. 
[0085] Geometric measures are constructed to estimate and bound the amplitude wanderings to facilitate detection. 
In particular we seek a measure for each mutation-recognizer centered (MRC-) hixel that is invariant to local degra- 
dation. Measure which can be expressed by multiple integrals of the form 

50 
55 

where Z denotes the set of mutations of interest. In other words, we determine the function f(z) under the condition 
that m(z) should be invariant with respect to all dispersions Also, up to a constant factor, this measure is the only 
one which is invariant under a group of motions in a plane. In principle, we derive deterministic analytical transformations 
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on each MRC-hixel., that map error-elliptic dispersion bound defined on«H 2 (the two dimension Euclidean space - i.e., 
oxel layout) onto measures defined ong?. The dispersion bound is given by 



[0086] Such a representation of uniqueness facilitates the rapid decimation of the search space. It is implemented 
using a filter constructed using measure-theoretic arguments. The transformation under consideration has its theoret- 
ical basis in the Palm Distribution Theory for point processes in Euclidean spaces, as well as in a new treatment in the 
problem of probabilistic description of MRC-hixel dispersion generated by a geometrical processes. Latter is reduced 
to a calculation of intensities of point processes. Recall that a point process in some product space E X F is a collection 
of random realizations of that space represented as {(ej, fj), | e } G. E, f } e F }. 

[0087] The Palm distribution, n of a translation (T n ) invariant, finite intensity point process in 91" is defined to the 
conditional distribution of the process. Its importance is rooted in the fact that it provides a complete probabilistic 
description of a geometrical process. 

[0088] In the general form, the Palm distribution can be expressed in terms of a Lebesgue factorization of the form 



[0089] Where n and A completely and uniquely determine the source distribution P of the translation invariant point 
process. Also, E P N* denotes the first moment measure of the point process and I_m is a probability measure. 
[0090] Thus a determination of n and A is needed which can uniquely encode the dispersion and amplitude wandering 
associated with the MRC-hixel. This is achieved by solving a set of equations involving Palm Distribution for each 
hybridization (i.e., mutation of interest). Each hybridization is treated as a manifestation of a stochastic point process 
in 9? 1 

[0091] In order to determine n and A we have implemented the following measure- theoretic filter: 
Determination of A 

using integral formulae constructed using the marginal density functions for the point spread associated with 
MRC-hixel(i f j) 

The oligonucleotide density per oxel p m (/j). p CR amplification protocol (o m ), fluorescence binding efficiency 
(T| m ) and imaging performance (o> m ) provide the continuous probability density function for amplitude wandering 
in the m-th MRC-hixel of interest. Let this distribution be given by (p m (/j), o m , r\ m , co m ). 

The method requires a preset binding dispersion limit to be provided to compute A. The second moment to 



Log^^l 2 ). 



E P N* = A L N X n 



the function 



(Pm(ij)- °m' ^rrr w /n> 



at SNR = 0 condition is used to provide the bound. 



Determination of fl 

Obtained by solving the inverse problem 



n = e*p 



where 




where x 1 and t 2 represent the normalized hybridization dispersion limits. 
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[0092] The number are empirically plugged in. The values of 0.1 and 0.7 are appropriate for, respectively, signifying 
loss of 10% - 70% hybridization. Also, 0 denotes the distribution of known point process. The form 1/(1+exp(p (...))) is 
employed herein to represent 0. 

[0093] The biomicroarray dispersion coefficients, hybridization process variability values and empirical variance are 
5 then converted to parameters at step 304 for use in NIF as follows. 
[0094] The NIF is represented by: 

Predicted State = f(current state, observation model, information uncertainty, information model) 

10 

[0095] Detailed equations are given below. 

[0096] In the biomicroarray context, NIF is an information-theoretic filter that predicts and estimates information about 
nonlinear state parameters (quality of observable) given nonlinear observations (e.g., post hybridization imaging) and 
nonlinear system dynamics (spatio-temporal hybridization degradation). The NIF is expressed in terms of measures 
15 of the amount of information about the observable (i.e., parameter of interest) instead of tracking the states themselves, 
it has been defined as the inverse covariance of the Kalman filter, where the information is in the Fisher sense, i.e, a 
measure of the amount of information about o, present in the observations Z k where the Fisher information matrix is 
the covariance of the score function. 

[0097] In a classical sense the biomicroarray output samples can be described by the nonlinear discrete-time state 
20 transition equation of the form: 

VDC(k)= f(VDC(k - -\)Mk ' + v(k) 

25 where VDC(k-1 ) is the state at time instant (k-1 ), 

0(k-1 ) is the input vector (embodied by dosage and/or therapy) 

v(k) is some additive noise; corresponds to the biomicroarray dispersion as computed by the Palm Generator 
functions above. 

VDC(k) is the state at time k, 

30 f(k,.,) is the nonlinear state transition function mapping previous state and current input to the current state. In 

this case it is the fractal mapping that provides the VDC coordinate at time k. 

[0098] The observations of the state of the system are made according to a non-linear observation equation of the 
form 

35 

z(k) = h(VDC(k)) + w(k) 

where z(k) is the observation made at time k 
VDC(k) is the state at time k, 
40 w(k) is some additive observation noise 

and h(.,k) is the current non-linear observation model mapping current state to observations, i.e., sequence-by- 
hybridization made at time k, 

v(k) and w(k) are temporally uncorrected and zero-mean. This is true for the biomicroarray in how protocol 
uncertainties, binding dynamics and hybridization degradation are unrelated and additive. The process and observation 
45 noises are uncorrected. 

E[v(i) w T (j)] = 0, V i,j. 

50 [0099] The dispersion coefficients together define the nonlinear observation model. 
[0100] The nonlinear information prediction equation is given by 

y(k\k - 1) = Y(k\k - 1)f(fc, VDC(k - 1|k - 1),u(/c - 1)) 



Y (k | k - 1) = [v/ t V// {k) + Q(k)]- 1 
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[0101] The nonlinear estimation equations are given by 

y(*l*)= y(«*-i) + '(*) 



Y(k\k)= Y(k\k - A) + \(k) 

where 

10 

l{k) = Vhl(k) f?" 1 (fr)V/> x (k) 

T -1 A 

15 /(*) = V/j^/c) R (/c)[v(/c) + Vh x {k) VDC{k\k - 1)] 

where 

20 v(k) = z(k)-h(x(k\k-1)) 

[0102] In this method NIF helps to bind the variability in the VDC coordinate mapping across sample to sample so 
that dosage and therapy effectiveness can be accurately tracked. 

[0103] The NIF is then applied to enhanced, fractal-filtered dot spectrogram at step 306 as follows. States being 
25 tracked correspond to post-hybridization dot spectrogram in this method. NIF computation as described above specifies 
the order interval estimate associated with a VDC point. It will explain and bound the variability in Viral load estimations 
for the same patient from laboratory to laboratory. 

[0104] The NIF also specifies how accurate each VDC coordinate is given the observation model and nucleotide set 
being analyzed. 

30 

Convergence Testing 

[0105] Referring again to FIG. 1, once any uncertainty is compensated, the VDC coordinates are renormalized at 
step 114. The renormalized VDC coordinates are patient specific and therapy specific. Alternately the coordinates 

35 could be virus / nucleotide marker specific. The NIF-compensated VDC coordinates are renormalized to the first diag- 
nostic sample point obtained using the biomicroan-ay. Thus a patient can be referenced to any point on the VDC. 
[0106] This renormalization step ensures that VDC properties are maintained, notwithstanding information uncer- 
tainties as indicated by the NIF correction terms. The approach is drawn from "renormalization-group" approach used 
for dealing with problems with many scales. In general the purpose of renormalization is to eliminate an energy scale, 

40 length scale or any other term that could produce an effective interaction with arbitrary coupling constants. The strategy 
is to tackle the problem in steps, one step for every length scale. In this method the renormalization methodology is 
abstracted and applied during a posteriori regularization to incorporate information uncertainty and sample-to-sample 
variations. 

[01 07] This is in contradistinction to current viral load measurement calibration methods that either generate samples 
45 with same protocol and same assumptions of uncertainty or use some constant correction term. Both existing ap- 
proaches skew the viral load readout so that measurements are actually accurate only in a limited -information" and 
"observability" context. This explains the large variations in readings from different laboratories and technicians for the 
same patient sample. 

[0108] Specifically, we include the dynamic NIF correction function to the gradient of the VDC at the sample point 
so normalized in a manner such that when the information uncertainty is null, the correction term vanishes. As discussed 
in the above steps, the NIF correction terms is actually derived from the noise statistics of the microarray sample. 

< VDC'(k,Q) > = VDC(k,e) + [V/v7F( V,/)J 

55 

where VNIF(Y,l) k denotes the gradient of nonlinear information prediction function. Under perfect observation 
model this term vanishes. 

[0109] Once initialized, the VDC coordinates are then updated at step 1 16 applying the IFS filter W[ ] on k+1th sample, 
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by 

VDC(k+1, 0)«- W[Biomicroarray Output, K+1); 

[0110J A direction convergence test is next performed at step 118 to determine whether the selected therapy has 
been effective. If convergence establishes that the viral load for the patient is moving in a direction representative of 
a lower viral load, then the therapy is deemed effective. The system is deemed to be converging toward a lower viral 
load if and only if: 



VDC(t k )-VDC(t k _,) l > a XVDC^-VDClu) 
VDCQ^-VDCit^l > A IVDC^-VDCQ^) 



[011 1] The above relationships needs to be monotonically persistent for at least two combinations of k and j. 
[0112] Also, datefk] - date[j] < k * characteristic time, / (in days) 
[01 1 3] Where k captures the population variability. Typically, k < 1 .2 . 

[0114] The peak VDC value is determined based on the VDC. The peak amplitude is an artifact of the specific pa- 
rameterization to the Fokker-Planck equation used in deriving the VDC. It is almost always derived independent of the 
specific sample. 

[0115] In connection with step 1 1 8, a VDC Shift factor A may be specified at which a dosage effectiveness decision 
and/or disease progression decision can be made. The VDC shift factor is applied to estimate the VDC curvature 
traversed between two measurements. 

[0116] If the system is deemed to be converging toward a lower viral load, an output signal is generated indicating 
that the therapy of interest is effective at step 120. If not, then the execution proceeds to step 122 wherein VDC scale 
matching is performed. A key assumption underlying this method is that movement along VDC is significant if and only 
if the sample points are with a constant multiple of temporal scale characterizing the VDC. This does not in any way 
preclude the pharmacological relevance associated with the datapoints. But complete pharmacological interpretation 
ofthe sample points is outside the scope ofthis method. The process is assumed to be cyclostationary or at a large 
time scale and two or more sample points have been mapped to VDC coordinates. The coordinates are then plugged 
into an analytic to estimate the empirical cycle time ( / ). This is implemented as described in the following sections. 
[0117] Again the empirical cycle time ( / ) is used to establish decision convergence. 

Scale Matching 

[0118] Select a forcing function of the form: 

= k. #? m cos cot 

where k is a constant and m is a small odd integer (m<7). 
[01 19] The phase space for this dynamical system is represented by: 




where 



< 1 fork>2 and j>0 
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erf m(x) = 



-1 ifx<-N 
erf(x) if|x|<SN 
1 x>N 



and erf m(.) denotes the error function. 

K is set to 1 ; where o< a, y, co < I are refer to constants. 

Let t emp denote the cycle-time-scale for this empirical system. 

If log e (T emp /T) > 1 (in step 10) then we claim that time-scales do not match. 

Time Scale Testing 

[0120] Next a determination has been made as to whether an effectiveness timescale has been exceeded at step 
124 by: 

Checking if a time step between successive sampling has exceeded T by determining if 

Time k+1 - Time k > T 

such that VDC(k+1 , 0)- VDC(k, 0) < £ where 
is set to 0.0001 and wherein 
T is given by 



^ 1 Tt 2*073,1/3) aJSl 
r=-an*os^I m —J 

B(1/3,1/3) represents the Beta function around the coordinates (1/3,1/3), We can actually use all B(1/2+1 ,1/2i+1) 
for i>1 and i<7. 

If Timek+i - Time k > T then output signal at step 126 indicating that either 
no change in viral load concluded, OR 
therapy deemed ineffective, OR 
dosage deemed suboptimal. 

[0121] If Time k+1 - Tlme k < T then process another sample by repeating all steps beginning with Step 4 wherein a 
dot spectrogram is generated for a new sample. 

[0122] Ifthe effectiveness time scale has been exceeded then a signal is output indicating that no determination can 
be as to whether the therapy of interest is effective. If the time scale is not exceeded, then execution returns no step 
106 for processing another sample. If available, and the processing steps are repeated. 

Alternative Implementations 

[0123] Details regarding a related implementation may be found in U.S. Patent 6,142,681 entitled "Method and Ap- 
paratus for Interpreting Hybridized Bioelectronic DNA Microarray Patterns Using Self Scaling Convergent Reverberant 
Dynamics". 

[01 24] The exemplary embodiments have been primarily described with reference to flow charts illustrating pertinent 
features of the embodiments. Each method step also represents a hardware or software component for performing 
the corresponding step. These components are also referred to software component for performing the corresponding 
step. These components are also referred to herein as a "means for" performing the step. It should be appreciated 
that not all components of a complete implementation of a practical system are necessarily illustrated or described in 
detail. Rather, only those components necessary for a thorough understanding of the invention have been illustrated 
and described in detail. Actual implementations may contain more components or, depending upon the implementation, 
may contain fewer components. 
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Claims 



1 A method for determining the effectiveness of a therapy by analyzing microarray output patterns generated from 
biological samples taken at different sampling times from a patient undergoing the therapy, said method compns.ng 
the steps of: 

generating a viral diffusion curve associated with a therapy of interest; 

mapping each of the output patterns representative of hybridization activity to respective coordinates on the 
viral diffusion curve using fractal filtering; 

determining a degree of convergence between the mapped coordinates on the viral diffusion curves; and 
determining whether the therapy of interest has been effective based upon the degree of convergence from 
one sample to another. 

2. The method of daim 1 wherein the step of generating the viral diffusion curve includes the steps of 
inputting parameters representative of viral load studies for the therapy of interest; 
generating a preliminary viral diffusion curve based upon the viral load studies; and 
calibrating a degree of directional causality in the preliminary viral diffusion curve to yield the viral diffusion 



curve. 



3 The method of claim 2 wherein the parameters representative of the viral load studies include one or more of 
baseline viral load (BVL) set point measurements at which detection is achieved, BVL at which therapy is recom- 
mended and viral load markers at which dosage therapy is recommended. 

4. The method of claim 2 wherein the step of generating a preliminary viral diffusion curve is performed by: 
selecting an equation representative of the viral diffusion curve; 

determining expectation and mean response parameters for use in parameterizing the equation selected to 
represent the vira! diffusion curve; and . 
parameterizing the equation selected to represent the viral diffusion curve to yield the preliminary viral diffusion 

curve. 

5 The method of claim 4 wherein the step of determining expectation and mean response parameters for use in 
parameterizing the equation selected to represent the viral diffusion curve comprises the steps of: 

performing PCR amplification; 

obtaining calibrated viral counts from the PCR amplification; 

determining the enhanced and normalized hybridization amplitude mean and variance values corresponding 
to the calibrated viral counts; and 

matching the enhanced and normalized hybridization amplitude mean and variance values to yield the expec- 
tation and mean response parameters. 

6. The method of claim 2 wherein the step of calibrating a degree of directional causality in the preliminary viral 
diffusion curve to yield the viral diffusion curve employs nonlinear information filtering. 

7. The method of claim 1 wherein the step of mapping each output pattern to the viral diffusion curve using fractal 
filtering comprises the steps of: 

generating a partitioned iterated fractal system IFS model representative of the output pattern; 

determining affine parameters for IFS model; and 

mapping the output pattern onto the viral diffusion curve using the IFS. 

8. The method of claim 1 wherein the step of determining whether the therapy of interest has been effective includes 
the steps of: 

evaluating the degree of convergence against a predetermined threshold of convergence; 

if the degree of convergence does not exceed the predetermined threshold, outputting a signal indicating that 

the therapy of interest has been effective; and 

if the degree of convergence exceeds the predetermined threshold, indicating a lack of effectiveness of the 
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therapy, then performing viral diffusion curve scale matching on the filtered vira! diffusion curve and then eval- 
uating whether a predetermined effectiveness time scale has been exceeded, and 

if the time scale has been exceeded, outputting a signal indicating the therapy has not been effective, and 
if the time scale has not been exceeded, then obtaining a new sample and repeating all steps using the pre- 
5 viously generated viral diffusion curve. 

9. The method of claim I further including the step of compensating for uncertainty in the filtered viral diffusion curves. 

10. The method of claim 9 wherein the step of compensating for uncertainty in the filtered viral diffusion curves is 
10 performed using non-linear information filtering. 

11. The method of claim 1 further including the step, performed before the mapping step, of enhancing the output 
pattern interferometrically. 

15 1 2. The method of claim 11 wherein the step of enhancing the output pattern interferometrically comprises the step of 
performing a nucleic acid amplification. 

13. The method of claim I wherein the output pattern is a quantized output pattern. 

20 14. The method of claim 13 wherein the quantized output pattern is a dot spectrogram. 

15. The method of claim 1 wherein the biological sample is selected from a group consisting of a DNA, RNA, protein, ' 
peptide-nucleic acid (PNA) and targeted nucleic amplification (TNA) samples. 

25 16. A system for determining the effectiveness of a therapy by analyzing microarray output patterns generated from 
biological samples taken at different sampling times from a patient undergoing the therapy, said system comprising: 

a viral diffusion curve generator for generating a viral diffusion curve associated with a therapy of interest; 
a mapping unit for mapping each of the output patterns representative of hybridization activity to respective 
30 coordinates on the viral diffusion curve using fractal filtering; 

convergence determination unit for determining a degree of convergence between the mapped coordinates 
on the viral diffusion curves for determining whether the therapy of interest has been effective based upon the 
degree of convergence from one sample to another. 

35 17. The system of claim 16 wherein the viral diffusion curve generator includes 

an input unit for inputting parameters representative of viral load studies for the therapy of interest; 
a preliminary viral diffusion curve generator for generating a preliminary viral diffusion curve based upon the 
viral load studies; and 

a calibration unit for calibrating a degree of directional causality in the preliminary viral diffusion curve to yield 
40 the viral diffusion curve. 

18. The system of claim 16 wherein the parameters representative of the viral load studies include one or more of 
baseline viral load (BVL) set point measurements at which detection is achieved, BVL at which therapy is recom- 
mended and viral load markers at which dosage therapy is recommended. 

45 

19. The method of claim 17 wherein the preliminary viral diffusion curve generator operates by 

selecting an equation representative of the viral diffusion curve; 

determining expectation and mean response parameters for use in parameterizing the equation selected to 
represent the viral diffusion curve; and 
50 parameterizing the equation selected to represent the viral diffusion curve to yield the preliminary viral diffu- 

sion curve. 

20. The system of claim 19 wherein the preliminary viral diffusion curve generator operates to determine expectation 
and mean response parameters for use in parameterizing the equation selected to represent the viral diffusion 

55 curve by: 

performing PCR amplification; 

obtaining calibrated viral counts from the PCR amplification; 
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determining the enhanced and normalized hybridization amplitude mean and variance values corresponding 
to the calibrated viral counts; and 

matching the enhanced and normalized hybridization amplitude mean and variance values to yield the expec- 
tation and mean response parameters. 

21. The system of claim 17 wherein the calibration unit employs nonlinear information filtering. 

22. The system of claim 16 wherein the mapping unit operates by: 

generating a partitioned iterated fractal system IFS model representative of the output pattern; 

determining affine parameters for IFS model; and 

mapping the output pattern onto the viral diffusion curve using the IFS. 

23. The system of claim 16 wherein the convergence determination unit operates by: 

evaluating the degree of convergence against a predetermined threshold of convergence; 

if the degree of convergence does not exceed the predetermined threshold, outputting a signal indicating that 

the therapy of interest has been effective; and 

if the degree of convergence exceeds the predetermined threshold, for indicating a lack of effectiveness of 
the therapy and for performing viral diffusion curve scale matching on the filtered viral diffusion curve and then 
evaluating whether a predetermined effectiveness time scale has been exceeded, and 
if the time scale has been exceeded, outputting a signal indicating the therapy has not been effective, and 
if the time scale has not been exceeded, obtaining a new sample and repeating all steps using the previously 
generated viral diffusion curve. 

24. The system of claim 16 further including a compensation unit for compensating for uncertainty in the filtered viral 
diffusion curves. 

25. The system of claim 24 wherein the compensation unit employs non-linear information filtering. 

26. The system of claim 16 further including an interferometric enhancement unit for enhancing the output pattern 
interferometrically. 

27. The system of claim 26 wherein the interferometric enhancement unit operates to perform a nucleic acid amplifi- 
cation. 

28. The system of claim 26 wherein the output pattern is a quantized output pattern. 

29. The system of claim 28 wherein the quantized output pattern is a dot spectrogram. 

30. The system of daim 26 wherein the biological sample is selected from a group consisting of a DNA, RNA, protein, 
peptide-nucleic acid (PNA) and targeted nucleic amplification (TNA) samples. 

31. The method according to one of the claims 1 to 15 wherein output patterns of a biochip are used as the microarray 
output patterns. 

32. The system according to one of the claims 16 to 30 for analyzing output patterns of a biochip as the microarray 
output patterns. 

Patentanspruche 

1. Verfahren zum Ermitteln der Wirksamkeit einer Therapie durch Analyse von aus zu unterschiedlichen Probeent- 
nahmezeitpunkten von einem der Therapie unterzogenen Patienten entnommenen biologischen Proben erzeugten 
Mikroarrayausgabemustem, wobei das Verfahren die Schritte umfasst: 

- Erzeugen einer einer im Interesse stehenden Therapie zugeordneten Virusdiffusionskurve, 

- Abbilden jedes der Hybridisierungsaktivitat angebenden Ausgabemuster auf entsprechende Koordinaten auf 
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der Virusdiffusionskurve unter Verwendung fraktaler Filterung, 

Ermitteln eines Konvergenzgrads zwischen den abgebildeten Koordinaten auf den Virusdiffusionskurven, und 
Ermitteln, ob die im Interesse stehende Tberapie wirksam war, basierend auf dem Konvergenzgrad von einer 
Probe zu einer Anderen. 

Verfahren nach Anspmch 1, bei dem der Schritt, die Virusdiffusionskurve zu erzeugen, die Schritte aufweist: 

- Eingeben von Parameter^ die Virusbelastungsstudien fur die im Interesse stehende Therapie angeben, 
Erzeugen einer voriaufigen Virusdiffusionskurve basierend auf den Virusbelastungsstudien, und 
Kalibrieren eines richtungsabhangigen Kausalitatsgrads in der voriaufigen Virusdiffusionskurve, urn die Virus- 
diffusionskurve zu erhalten. 

Verfahren nach Anspruch 2, bei dem die die Virusbelastungsstudien angebenden Parameter von Basislinienvirus- 
belastungs-(BVL)-Sollwertmessungen, bei denen Detektion erreicht wird, BVL, bei der Therapie empfohlen wird, 
und Wusbelastungsmarkern, bei denen Dosistherapie empfohlen wird, einen oder mehr aufweisen. 

Verfahren nach Anspruch 2, bei dem der Schritt, eine vorlaufige Virusdiffusionskurve zu erzeugen, durchgefuhrt 
wird durch: 

- Wahlen einer die Virusdiffusionskurve angebenden Gleichung, 

Ermitteln von Erwartungs- und mittleren Antwortparametem zur Verwendung beim Parametrisieren der zur 
Angabe der Virusdiffusionskurve gewahlten Gleichung und 

Parametrisieren der zur Angabe der Virusdiffusionskurve gewahlten Gleichung, urn die vorlaufige Virusdiffu- 
sionskurve zu erhalten. 

Verfahren nach Anspmch 4, bei dem der Schritt, Erwartungs- und mittlere Antwortparameter zur Verwendung beim 
Parametrisieren der zur Angabe der Virusdiffusionskurve gewahlten Gleichung zu ermitteln, die Schritte umfasst: 

Durchfuhren einer PCR-Verstarkung, 

Erhalten von kalibrierten Viruszahlungen aus der PCR-Verstarkung, 

Ermitteln der den kalibrierten Viruszahlungen entsprechenden Mittel- und Varianzwerten der verstarkten und 
normalisierten Hybridisierungsamplituden, und 

- Abgleichen der Mittel- und Varianzwerte der verstarkten und normalisierten Hybridisierungsamplituden, urn 
die Erwartungs- und mittleren Antwortparameter zu erhalten. 

Verfahren nach Anspruch 2, bei dem der Schritt, einen richtungsabhangigen Kausalitatsgrad in der voriaufigen 
Virusdiffusionskurve zu kalibrieren, urn die Virusdiffusionskurve zu erhalten, nicht lineare Informationsfilterung 
verwendet. 

Verfahren nach Anspruch 1 , bei dem der Schritt, jedes Ausgabemuster auf die Virusdiffusionskurve unter Verwen- 
dung fraktaler Filterung abzubilden, die Schritte umfasst: 

Erzeugen eines das Ausgabemuster angebenden, partitionierten iterierten Fraktalsystem-IFS-Modell, 
Ermitteln affiner Parameter fur das IFS-Modell, und 

Abbilden des Ausgabemusters auf die Virusdiffusionskurve unter Verwendung des IFS. 

Verfahren nach Anspruch 1 , bei dem der Schritt, zu Ermitteln, ob die im Interesse stehende Therapie wirksam war, 
die Schritte aufweist: 

Evaluieren des Konvergenzgrads gegenuber einem vorbestimmten Konvergenzgrenzwert, 

- Ausgeben eines Signals, das angibt, dass die im Interesse stehende Therapie wirksam war, wenn der Kon- 
vergenzgrad den vorbestimmten Grenzwert nicht uberschreitet, und 

Angeben eines Fehlens der Wirksamkeit der Therapie, wenn der Konvergenzgrad den vorbestimmten Grenz- 
wert uberschreitet, wobei dann eine Skalierungsanpassung der Virusdiffusionskurve auf die gefllterte Virus- 
diffusionskurve durchgefuhrt und dann evaluiert wird, ob ein vorbestimmter WirksamkeitszeitmafJstab uber- 
schritten wurde, und 

Ausgeben eines Signals, das angibt, dass die Therapie nicht wirksam war, wenn der Zeitmafistab uberschritten 
wurde, und 
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- Erhalten einer neuen Probe und Wiederholen aller Schritte unter Verwendung der zuvor erzeugten Virusdif- 
fusionskurve dann, wenn der ZeitmaBstab nicht uberschritten wurde. 

9. Verfahren nach Anspruch 1 , ferner den Schritt aufweisend, Unscharfen in den gefilterten Virusdiffusionskurven zu 
5 kompensieren. 

10. Verfahren nach Anspnjch 9, be! dem der Schritt, Unscharfen in den gefilterten Virusdiffusionskurven zu kompen- 
sieren, unter Verwendung nicht linearer Informationsfilterung durchgefuhrt wird. 

w 11. Verfahren nach Anspruch 1 , femer den vor dem Abbildungsschritt durchgefuhrten Schritt aufweisend, das Ausga- 
bemuster interferometrisch zu verstarken. 

12. Verfahren nach Anspruch 11 , bei dem der Schritt, das Ausgabemuster interferometrisch zu verstarken, den Schritt 
umfasst, eine Nukleinsaureverstarkung durchzufuhren. 

15 

13. Verfahren nach Anspruch 1, bei dem das Ausgabemuster ein quantisiertes Ausgabemuster ist. 

14. Verfahren nach Anspruch 13, bei dem das quantisierte Ausgabemuster ein Punktspektrogramm ist. 

20 15. Verfahren nach Anspruch 1, bei dem die biologische Probe aus einer Gruppe ausgewahlt wird, die aus DNA-, 
RNA-, Protein-, Peptid-Nukleinsaure-(PNA)- und Target-Nukleidverstarkungs-(TNA)-Proben besteht. 

16. System zum Ermitteln der Wirksamkeit einer Therapie durch Analyse von aus zu unterschied lichen Probenent- 
nahmezeitpunkten von einem der Therapie unterzogenen Patienten erhaltenen biologischen Proben erzeugten 
25 Mirkoarrayausgabemustem, wobei das System umfasst: 

- einen Virusdiffusionskurvengenerator zum Erzeugen einer einer im Interesse stehenden Therapie zugeord- 
neten Virusdiffusionskurve, 

- eine Abbildungseinheit zum Abbilden jedes der Hybridisierungsaktivitat angebenden Ausgabemuster auf ent- 
30 sprechende Koordinaten auf der Virusdiffusionskurve unter Verwendung fraktater Filterung, 

- einer Konvergenzermittlungseinheit zum Ermitteln eines Konvergenzgrads zwischen den abgebildeten Koor- 
dinaten auf den Virusdiffusionskurven, urn zu ermitteln, ob die im Interesse stehende Therapie wirksam war, 
basierend auf dem Konvergenzgrad von einer Probe zu einer Anderen. 

35 17. System nach Anspruch 16, bei dem der Virusdiffusionskurvengenerator aufweist: 

- eine Eingabeeinheit zum Eingeben von Parametem, die Virusbelastungsstudien uber die im Interesse ste- 
hende Therapie angeben, 

- einen Generator fur eine vorlaufige Virusdiffusionskurve, urn eine vorlaufige Virusdiffusionskurve basierend 
40 auf den Virusbelastungsstudien zu erzeugen, und 

- eine Kalibrierungseinheit zum Kalibrieren eines richtungsabhangigen Kausalitatsgrads in der vorlaufigen Vi- 
rusdiffusionskurve, urn die Virusdiffusionskurve zu erhalten. 

1 8. System nach Anspruch 1 6, bei dem die die Virusbelastungsstudien angebenden Parameter von Basislinienvirus- 
45 belastungs-(BVL)-Sollwertmessungen, bei denen Detektion erreicht wird, BVL, bei der Therapie empfohlen wird, 

und Virusbelastungsmarkem, bei dem eine Dosistherapie empfohlen wird, eine Oder mehr aufweisen. 

19. Verfahren nach Anspruch 17, bei dem der Generator fur eine voriaufige Virusdiffusionskurve arbeitet durch: 

50 - Wahlen einer die Virusdiffusionskurve angebenden Gleichung, 

- Ermitteln von Erwartungs- und mittleren Antwortparametem zur Verwendung beim Parametrisieren der zur 
Angabe der Virusdiffusionskurve gewahlten Gleichung, und 

- Parametrisieren der zur Angabe der Virusdiffusionskurve gewahlten Gleichung, urn die vorlaufige Virusdiffu- 
sionskurve zu erhalten. 



55 



20. System nach Anspruch 19, bei dem der Generator fur eine vorlaufige Virusdiffusionskurve zum Ermitteln von 
Erwartungs- und mittleren Antwortparametem zur Verwendung beim Parametrisieren der zur Angabe der Virus- 
diffusionskurve gewahlten Gleichung arbeitet durch: 
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Durchfuhren einer PCR-Verstarkung, 

Erhalten kalibrierter Vlruszahlungen aus der PCR-Verstarkung, 

Ermitteln der den kalibrierten Viruszahlungen entsprechenden Mittel- und Varianzwerten der verstarkten und 
normal isierten Hybridisierungsamplituden, und 

Abgleichen der Mittel- und Varianzwerte der verstarkten und normal isierten Hybridisierungsamplituden, urn 
die Erwartungs- und mittleren Antwortparameter zu erhalten. 

21. System nach Anspruch 17, bei die Kalibrierungseinheit nicht lineare Infonmationsfilterung verwendet. 

22. System nach Anspruch 16, bei dem die Abbildungseinheit arbeitet durch: 

Erzeugen eines das Ausgabemuster angebenden, partition ierten iterierten Fraktalsystem-IFS-Modells, 
Ermitteln affiner Parameter fur das IFS-Modell, und 

Abbilden des Ausgabemusters auf die Virusdiffusionskurve unter Verwendung des IFS. 

23. System nach Anspruch 16, bei dem die Konvergenzermittlungseinheit arbeitet durch: 

Evaluieren des Konvergenzgrades gegenuber einem vorbestimmten Konvergenzgrenzwert, 

Ausgeben eines Signals, das angibt, dass die im Interesse stehende Therapie wirksam war, wenn der Kon- 

vergenzgrad den vorbestimmten Grenzwert nicht uberschreitet, und 

Angeben eines Fehlens der Wirksamkeit der Therapie und Durchfuhren einer Skalierungsanpassung der Vi- 
rusdiffusionskurve an die gefilterte Virusdiffusionskurve und dann Evaluieren, ob ein vorbestimmter Wirksam- 
keitszeitmalistab uberschritten wurde, wenn der Konvergenzgrad den vorbestimmten Grenzwert uberschrei- 
tet, und 

Ausgeben eines Signals, das angibt, dass die Therapie nicht wirksam war, wenn der Zeitmafistab uberschritten 
wurde, und 

Erhalten einer neuen Probe und Wiederholen aller Schritte unter Verwendung der zuvor erzeugten Virusdif- 
fusionskurve, wenn der ZeitmafJstab nicht uberschritten wurde. 

24. System nach Anspruch 16, ferner mit einer Kompensationseinheit, urn Unscharfen in den gefilterten Virusdiffusi- 
onskurven zu kompensieren. 

25. System nach Anspruch 24, bei dem die Kompensationseinheit nicht lineare Informationsfilterung verwendet. 

26. System nach Anspruch 16, ferner mit einer interferometrischen Verstarkungseinheit, urn das Ausgabemuster in- 
terferometrisch zu verstarken. 

27. System nach Anspruch 26, bei dem die interferometrische Verstarkungseinheit arbeitet, um eine Nukleidsaure- 
verstarkung durchzufuhren. 

28. System nach Anspruch 26, bei dem das Ausgabemuster ein quantisiertes Ausgabemuster ist. 

29. System nach Anspruch 28, bei dem das quantisierte Ausgabemuster ein Punktspektrogramm ist. 

30. System nach Anspruch 26, bei dem die biologische Probe von einer Gruppe ausgewahlt ist, die aus DNA-, RNA-, 
Protein-, Peptid-Nukleidsaure-(PNA)- und Targetnucleidverstarkungs-(TNA)-Proben besteht. 

31. Verfahren nach einem der Anspruche 1 bis 15, bei dem Ausgabemuster eines Biochips als die Mikroarrayausga- 
bemuster verwendet werden. 

32. System nach einem der Anspruche 16 bis 30, um Ausgabemuster eines Biochips als die Mikroarrayausgabemuster 
zu analysieren. 

Revendications 

1. Precede de determination de I'efficacite d'une therapie par analyse de modeles de resultats de micro-depots ge- 
neres a partir d'echantillons biologiques preleves a differents points d'echantillonnage chez un patient recevant 
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cette therapie, ledit procede comprenant les etapes consistant a : 

generer une courbe de diffusion virale associee a la therapie concemee ; 

cartographier, sous forme de coordonnees respectives, chacun des modeles de resultats representatifs de 
I'activite d'hybridisation sur la courbe de diffusion virale a I'aide d'une technique de filtration fractale ; 
determiner un degre de convergence entre les coordonnees respectives des differents modeles de resultats 
cartographies sur les courbes de diffusion virale ; et 

determiner I'efficacite de la therapie concemee d'apres le degre de convergence d'un echantillon a un autre. 

Procede selon la revendication 1, dans lequel I'etape consistant a generer la courbe de diffusion virale inclut les 
etapes consistant a : 

entrer les parametres representatifs d'etudes de charge virale pour la therapie concemee ; 

generer une courbe de diffusion virale preliminaire sur la base des etudes de charge virale ; et 

calibrer un degre de causalite directionnelle a partir de la courbe de diffusion virale preliminaire pour donner 

la courbe de diffusion virale. 

Procede selon la revendication 2, dans lequel les parametres representatifs des etudes de charge virale incluent 
un ou plusieurs points de mesure predefinis de la charge virale de base (CVB) auxquels la detection est possible, 
CVB a laquelle la therapie est recommandee et marqueurs que constitue la charge virale aux quantites desquels 
la therapie est recommandee. 

Procede selon la revendication 2, dans lequel I'etape consistant a generer une courbe de diffusion virale est con- 
duce en : 

choisissant une equation representative de la courbe de diffusion virale ; 

determinant des parametres de valeur attendue et de reponse moyenne utilises dans le parametrage de liqua- 
tion choisie pour representer la courbe de diffusion virale ; et 

parametral I'equation choisie pour representer la courbe de diffusion virale afin de donner la courbe de dif- 
fusion virale preliminaire. 

Procede selon la revendication 4, dans lequel I'etape consistant a determiner des parametres de valeur attendue 
et de reponse moyenne utilises dans le parametrage de I'equation choisie pour representer la courbe de diffusion 
virale comprend les etapes consistant a : 

proceder a une amplification par PCR ; 

obtenir une quantite de virus calibree a partir de I'amplification par PCR ; 

determiner les valeurs de moyenne et de variance de I'hybridisation accrue et normalisee correspondant a la 
quantite de virus ; et 

harmoniser les valeurs de moyenne et de variance de I'hybridisation accrue et normalisee correspondant a 
la quantite de virus pour donner les parametres de valeur attendue et de reponse moyenne. 

Procede selon la revendication 2, dans lequel I'etape consistant a calibrer un degre de causalite directionnelle a 
partir de la courbe de diffusion virale preliminaire pour donner la courbe de diffusion virale utilise une technique 
de filtration de I'information non lineaire. 

Procede selon la revendication 1 , dans lequel I'etape consistant a cartographier chacun des modeles de resultats 
sur la courbe de diffusion virale a I'aide d'une technique de filtration fractale comprend les etapes consistant a : 

generer un modele de systeme fractal itere (IFS) partitionne representatif du modele de resultat ; 
determiner des parametres affins pour le modele IFS ; et 

cartographier le modele de resultat par dessus la courbe de diffusion virale a I'aide du modele IFS. 

Procede selon la revendication 1 , dans lequel I'etape consistant a determiner I'efficacite de la therapie concemee 
comprend les etapes consistant a : 

evaluer le degre de convergence par rapport a un seuil de convergence predetermine ; 

si le degre de convergence ne depasse pas le seuil predetermine, produire un signal indiquant que la therapie 
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convergence depasse le seuil predetermine, indiquer le manque d'efficacite de la therapie. puis 
hamSstl »a coume de diffusion virale d'apres la courbe de diffusion virale filtree e« defender 
si rechellede temps d'efficacite predeterminee est depassee;et 

si I'echelle de temps est depassee. produire un signal indiquant que la therapie n est pas efficace, e 

s! .■fchelle de temps n'est pas depassee, obtenir alors un nouve. echantillon et repeter toutes les etapes a 

I'aide de la courbe de diffusion virale pr6c6demment g§n6r6e. 

9. Precede selon la revendication 1 incluant en outre I'etape consistant a compenser I'incertitude associee aux cour- 
bes de diffusion virale filtrees. 

10. Precede selon .a revendication 9, dans lequel I'etape consistant a compenser incertitude associ^ , aux courbes 
de diffusion virale filtrees est conduite a I'aide d'une technique de filtration de I informat.on non l.nea,re. 

11. Precede selon la revendication 1 incluant en outre I'etape consistant a accrottre le modele de resultat interfero- 
rrtetriquement, conduite avant I'etape de cartographie. 

12. Precede selon la revendication 11 , dans lequel I'etape consistant a accrottre le modele de resultat interferometri- 
quement comprend I'Stape d'amplification d'acides nucl6iques. 

13. Precede selon la revendication 1 , dans lequel le modele de resultat est un modele de resultat quantifie. 

14. Precede selon la revendication 13, dans lequel le modele de resultat quantifie est un spectrogramme matriciel. 

15 Precede selon la revendication 1. dans lequel I'echantillon biologique est choisi parmi un groupe consistant en 
ADN ARN June pmteine, un acide peptide nucleique (APN) et des echantillons d'amplification nucle.que cb.ee 
(ANC). 

16 Svsteme de determination de I'efficacite d'une therapie par analyse de modules de resultats de micro-depots 

^a partt d'chantiHons biologiques preleves a differents points d'echantillonnage Chez un pat.en. recevant 
cette therapie. ledit systeme comprenant : 

un g6nerateur de courbe de diffusion virale permettant de generer une courbe de diffusion virale associee a 

unrSfdrcTrtograpnie permettant de cartographier. sous forme de coordonnees respectives ^un des 
modeles de resultats representatifs de I'activite d'hybridisatien sur la courbe de d,ffus,on v.rale a I a.de d une 

technique de filtration fractale ; . . 

une unite de determination de la convergence permettant de determiner un degre de ^^^S 
coordonnees respectives des differents modeles de resultats cartograph.es sur les courbes de ; dtfus on virale 
afin de determiner I'efficacite de la therapie concemee d'apres le degre de convergence d'un echant.llon a un 



autre. 



17 Svsteme selon la revendication 16. dans lequel le gen6rateur de courbe de diffusion virale inclut 

17. S y s, ^ n e e s 2 ,a ; e ntr6e pemlettant d . entrer les param6tres representatifs des etudes de charge vrale pour la 

SnSteurde courbe de diffusion virale preliminaire permettant de gender une courbe de diffusion virale 

preliminaire d'apr6s les etudes de charge virale ; et 

une unite de calibrage permettant de calibrer un degre de causalite direcfionnelle a part.r de la courbe de 
diffusion virale preliminaire pour donner la courbe de diffusion virale. 

18 Svsteme selon la revendication 16. dans lequel les parametres representatifs des etudes de charge virale incluent 
un ou plusieurs points de mesure pr^definis de la charge virale de base (CVB) auxquels la detect™ est poss-b^ 
CVB d , laquelle la therapie est recommand6e et marqueurs que constitue la charge v.rale aux quant.tes desquels 
la therapie est recommandee. 

19. Procede selon la revendication 17, dans lequel le generateur de courbe de diffusion virale preliminaire opere en : 
choisissant une Equation representative de la courbe de diffusion virale ; 
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determinant des parametres de valeur attendue et de reponse moyenne utilises dans le parametrage de I'equa- 
tion choisie pour representer la courbe de diffusion virale ; et 

parametral I'equation choisie pour representer la courbe de diffusion virale afin de donner la courbe de dif- 
fusion virale preliminaire. 

5 

20. Systeme selon la revendication 19, dans lequel le generateur de courbe de diffusion virale preliminaire determine 
les parametres de valeur attendue et de reponse moyenne utilises dans le parametrage de I'equation choisie pour 
representer la courbe de diffusion virale en : 

10 procedant a une amplification par PCR ; 

obtenant une quantite de virus calibree a partir de ('amplification par PCR ; 

determinant les valeurs de moyenne et de variance de I'hybridisation normalisee et accrue correspondant a 
la quantite de virus ; et 

harmonisant les valeurs de moyenne et de variance de I'hybridisation accrue et normalisee correspondant a 
15 la quantite de virus pour donner les parametres de valeur attendue et de reponse moyenne. 

21 . Systeme selon la revendication 1 7, dans lequel I'unite de calibrage utilise une technique de filtration de I'information 
non lineaire. 

20 22. Systeme selon la revendication 16, dans lequel I'unite de cartographie opere en : 

generant un modele de systeme fractal itere (IFS) partitionne representatif du modele de resultat ; 
determinant des parametres affins pour le modele IFS ; et 

cartographiant le modele de resultat par dessus la courbe de diffusion virale a I'aide 
25 du modele IFS. 

23. Systeme selon la revendication 16, dans le quel I'unite de determination de la convergence opere en : 

evaluant le degre de convergence par rapport a un seuil de convergence predetermine ; 
30 si le degre de convergence ne depasse pas le seuil predetermine, en produisant un signal indiquant que la 

therapie concemee est efficace ; et 

si le degre de convergence depasse le seuil predetermine, en indiquant un manque d'efficacite de la therapie, 
puis en harmonisant I'echelle de la courbe de diffusion virale d'apres la courbe de diffusion virale filtree et en 
determinant si I'echelle de temps d'efficacite predeterminee est depassee ; et 
35 si I'echelle de temps est depassee, en produisant un signal indiquant que la therapie concemee n'est pas 

efficace ; et 

si I'echelle de temps n'est pas depassee, en obtenant alors un nouvel echantillon et en repetant toutes les 
etapes a I'aide de la courbe de diffusion virale precedemment generee. 

40 24. Systeme selon la revendication 16 incluant en outre une unite de compensation de I'incertitude associee aux 
courbes de diffusion virale filtrees. 

25. Systeme selon la revendication 24, dans lequel I'unite de compensation utilise une technique de filtration de I'in- 
formation non lineaire. 

45 

26. Systeme selon la revendication 1 6 incluant en outre une unite d'accroissement interferometrique permettant d'ac- 
croitre le modele de resultat interferometriquement. 

27. Systeme selon la revendication 26, dans lequel le fonctionnement de I'unite d'accroissement interferometrique 
50 comprend I'amplification d'acides nucleiques. 

28. Systeme selon la revendication 26, dans lequel le modele de resultat est un modele de resultat quantifie. 

29. Systeme selon la revendication 28, dans lequel le modele de resultat quantifie est un spectrogramme matriciel. 

55 

30. Systeme selon la revendication 26, dans lequel I'echantillon biologique est choisi parmi un groupe consistent en 
ADN, ARN, une proteine, un acide peptide nucleique (APN) et des echantillons d'amplification nucteique ciblee 
(ANC). 
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31. Precede selon Tune quelconque des revendications 1 a 15, dans lequel les modules de resuitats d'une biopuce 
sont utilises comme model es de resuitats de micro-depots. 

32. Systeme d'analyse de modeles de resuitats d'une biopuce comme modeles de resuitats de micro-depots selon 
5 Tune quelconque des revendications 16 a 30. 
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FIGURE 5 
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